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ABSTRACT 

We report the results of our study of fractional entropy enhancement in the intracluster medium 
(ICM) of the clusters from the representative XMM-Newton cluster structure survey (REXCESS). We 
compare the observed entropy profile of these clusters with that expected for the ICM without any 
feedback, as well as with the introduction of preheating and cooling. We make the first estimate of the 
total, as well as radial, non-gravitational energy deposition upto 7-500 for this large, nearly fiux-limited, 
sample of clusters. We find that the total energy deposition corresponding to the entropy enhancement 
is proportional to the cluster temperature (and hence cluster mass) . The energy deposition, scaled by 
Tsp, per particle as a function of gas mass shows a similar profile in all clusters, with its being more 
pronounced in the central region than in the outer region. Our results support models of entropy 
enhancement through AGN feedback. 
Subject headings: galaxies: clusters : general — X-rays: galaxies : clusters 



1. INTRODUCTION 

Models of structure formation in the universe have 
been successful in predicting the average properties of 
galaxy clusters. Using these characteristics, such as 
gas temperature, XRay luminosity, SZ-fiux and rich- 
ness, it is possible to draw cosmol ogical conclus ions from 
survey s of galaxy clust e rs (e. g. . IReiprich fc B ohringcn 
(l200a): IVi khhnine t al.l (1200 9^: 'Gladder s et aJ.i f|200% 
Khedekar. Maium dar fc Das (2010); Ro zo et al.l |2010l ): 



Sehgal et al.l (I2011D : iBensMTet all ([20111). The detailed 
properties of the intracluster medium (ICM) however, 
need more input to the physics of baryonic gas than its 
falling i nto a dark mat t er po te ntial (e.g. | Shaw et al 



([20l0() : iBattagUa etliD (I20TTI): iTrac. Bode fc Ostriker 
(|20Tll) : iChaudhuri fc Maiumdaii (|201lD '). It is be- 



lieved that feedback from galaxies, including active 
galactic nuclei (AGNs), and/or radiative cooling of 
the ICM gas, modify the X-ray prope rties of the gas 
(see iMcNamara fc NulsenI (|2007l l2012| )). These non- 
gravitational processes tend to increase the entropy of the 
ICM gas, thereby making it tenuous, and consequently, 
under-luminous in X-rays, especially in low temperature 
(and mass) clusters. 
Recent observations of profiles of entropy (defined as 

K = ksT /nj , where n^ is the electron number density 
and /cb, the Boltzmann constant Q) allow one to compare 
them with theoretically expected profiles with or without 
feedback, and allow one to determine the nature and de- 
gree of feedback. Entropy as defined above is well suited 
for this sort of analysis as it is a record of the accretion of 
gas into the cluster, as well as the modificati ons shaped 
by the processes of gas cooling and feedback. IVoit et aTl 
(|2005| ) had shown that in the absence of any feedback and 
cooling processes, simulations tend to predict a power- 
law radial profile for the entropy outside the core, with 
a scaling K oc r^'^. 
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^ We write K for the entropy popularly defined in X-ray litera- 
ture, and denote the thermodynamic entropy as S 



Since entropy per particle is a Lagrangian quantity, it 
is more sensible to study the distribution of entropy not 
with the radial distance, but with gas mass, taking into 
accoui it the move r nent o f gas shells due to a change in en- 
tropy. iVoit et ahl ()2005D suggested the comparison of en- 
tropy as a function of gas mass in order to determine the 
enhance ment of entropy fron i non-g ravitational processes 
(see also lNath fc Majumdarl ()2011[ )). Here, we study the 
entropy profiles of clusters from the REXCESS sample 
and compare with the baseline profiles of ICM w ithout 
any feedback. For these clusters. iPratt et al.l ()2010[) stud- 
ied the radial entropy profiles and after comparing with 
the initial profile, found that entropy enhancement is ev- 
ident in the inner radii, and that it extends up to a large 
radii for low mass systems, while large mass clusters do 
not show entropy deviation at very large radii. In this 
Letter, we focus on the entropy profiles as functions of 
gas mass contained inside a given shell. We determine 
TdK/K , where dK/K is the fractional deviation of the 
observed entropy from the benchmark theoretically cal- 
culated entropy, which is a measure of the energy depo- 
sition per unit gas particle, and investigate the profile 
of this energy deposition for low and high temperature 
clusters. 

We adopt a ACDM cosmology with Hq = 70 km s^^ 
Mpc-i, ^M = 0.3 and Oa = 0.7. 

2. THE CLUSTER SAMPLE 

The REXCESS survey (|Bohringer et all [2001 uses 
the REFLEX cluster catalog as a parent sample. RE- 
FLEX is a nearly complete fiux limited cluster sam- 
le, covering 4.24 ster in the southern extragalactic sky 
Bohringer et al.|[2004D . This sample consists of 31 lo- 
cal clusters in the redshift range z ^ 0.2. The clus- 
ters arc selected on the basis of their X-Ray luminos- 
ity, Lx = (0.407-20) X 1044 h^^ erg s~^ in the 0.1-2.4 
keV band, with a homogeneous coverage in the chosen 
luminosity range, and no preference for any morphol- 
ogy type. The selected luminosity range provides clus- 
ters with a temperature ^ 2 kcV, and does not include 
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galaxy groups. As iPratt et alj (l201Clf ) have noted, the 
properties of the REXCESS sample allow one to study 
the variation of entropy profiles across a range of clus- 
ter masses, especially because the distances were chosen 
such that rsoo fell within the XMM-Newton field of view, 
which increased the precision of measurements at large 
radii. They also subdivided the sample into cool-core and 
non cool-core systems, defining the clusters with central 
density E{z)~''^ne.o > 4 x 10~^ cm~^ as cool-core systems 
{E{z) being the ratio of the Hubble constant at redshift 
z to its present value). 

In this work, we use the entropy profiles of 25 clus- 
ters from the who le REXCESS sample of 31 clusters (see 
IPratt et all ()2010D , their Tablel) . We use only those clus- 
ters with data at a minimum of 5 radial points outside 
the core radii, thus excluding clusters number 2, 13, 23, 
25 & 27 (ordered top-to-down respectively in the table). 
We also leave out cluster number 14 whose errors on ob- 
served entropy far exceed those of other clusters. 

3. ENTROPY PROFILES 

3.1. Initial entropy - radial profile 

In order to assess the entropy enhancement in observed 
clusters, we first discuss the pr ofile expected wi thout 
any non-gravitational processes. IVoit et all (|2005| ) pre- 
sented an analytic form for the baseline entropy profile 
which they obtained by analyzing the entropy profiles of 
clusters from non-radiative simulations. Their simulated 
SPH profiles, when fitted in the 0.1-1 r2oo range, scat- 
ter about a median scaled profile described by a baseline 
power law relation. 



K{r) 
K20Q 
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'^200 



1.1 
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with approximately 20% dispersion. These profiles how- 
ever are found to flatten within a radius of i? < 0.2 
rsoo, and they find that the agreement of the above fit 
with both their SPH and AMR simulations is better than 
^ 10% beyond a radius of 0.1 r2oo- 

3.2. Initial entropy profile with gas mass 

In this paper, we would like to study the entropy as 
a function of gas mass. In order to calculate the ini- 
tial entropy profile as a function of gas mass, we use 
the initial radial entropy profile, in conjunction with the 
assumption of hydrostatic equilibrium. We assume the 
Nava rro-Frenk- White (NEW) profile for the dark matter 
halo (jNavarro. Frenk fc White! |1997[ ) . For the concen- 
tration parameter C500 — 3.2 that we adopip, the corre- 



spond ing relation for eqn [T]at r5oo becomes, (jPratt et al.l 
[20T0h 

|W =1.42(^)11 (2) 

.f^SOO '"500 

The equation of hydrostatic equilibrium can be written 
as, 



dr 



-Pg- 



GM{< r) 



Pn 
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^ This v alue is measured for a mor phologically rel axed cluster 
sampl e by IPointecouteau et ai.l I I2005I ). also used by IPratt et al.l 
lf20Tg'l . 




Fig. 1. — Upper Panel - This plot shows the ratio of i^'th/i^soo 
as a function of gas mass Mg for all the clusters. Green lines refer 
to the lowest temperature clusters {Tsp < 3.5 keV), red for the 
intermediate temperature clusters (3.5 < Tsp < 5 keV), and blue 
lines are for the largest temperature clusters {Tsp > 5 keV) . Lower 
Panel - Kabs/K^oo is plotted against Mg/M^oo- Color scheme is 
the same as above. Dotted lines show observations below 0.1r200- 
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where Pg = iigkBT, is the gas pressure. For boundary 
conditions, we set the total gas fraction inside r^r to be 
the universal baryon fraction, fg = flb/flm- Eqns [2] & 
[3] are solved for the pressure profile Pg] this gives the 
density profile and hence Mg. One can then invert Mg 
to get K{Mg). 

We show the results in Figure [T] (Upper Panel) that 
plots the theoretical profiles of Kth scaled by /ism , the 
characteristic entropy (eq 3 in IPratt et al.l ()2010D ). for 
clusters of different temperature bins. The figure shows 
that Kth{Mg)/K5oo is self similar to a good approxi- 
mation. We have fitted the profile with a simple pa- 
rameterization, Kth{Mg)/K5oo = A{Mg/M5oo)", in the 
range 0.1r2oo — ?'500- The scatters in the slopes and nor- 
malization for different values of M500 are found to be 
small. The index a = 0.81 ± 0.05 for the whole sam- 
ple. The value of the parameter A was determined to be 
6.09 ±0.86. 
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TABLE 1 

Mean values of parameters in the range 0.1r2oo — ^500 
(excluding core): observed entropy-gas mass relation 



Sample 


Ao 


Bo 


do 


Total sample (25 clusters) 


0.23 ±0.43 


9.59 ± 7.54 


0.67 ± 0.47 


eool-core (9) 


0.14 ±0.61 


9.07 ±6.18 


0.63 ±0.61 


Non eool-core 


0.28 ±0.30 


9.89 ±8.38 


0.69 ± 0.40 



3.3. Observed entropy - radial profile 

IPratt et all (I20T0I ) have fit the REXCESS data to the 
form, 



/f(r)==A-o + /vioo[r/100kpc]' 



(4) 



where Kq is interpreted as the excess of core entropy 
above the best fitting power law at large radii. They 
scaled the quantities to rsoo , the effective limiting radius 
for high quality observations from XMM-Newton and 
Chandra, which they estimated iteratively from the up- 
dated calibration of the M500 — Yx relation, by including 
the REXCESS data for morphologically relaxed systems. 
Interior to rsoo, the observed entropy is always higher 
than the baseline prediction. At rsoo, they find the me- 
dian dimensionless entropy is i^(r5oo)/^5oo = 1-70 ± .35 
and that this is higher than but consistent with the base- 
line prediction. 

3.4. Observed entropy profiles with gas mass 

Next, we express the observed entropy profiles 
Kobs/K^oQ of REXCESS clusters as a function of 
Mg/Msoo- Figure [1] (Lower Panel) shows these profiles 
for all 25 clusters in our sample in different temperature 
ranges respectively. We fit the profiles by an expression 
of the form Kobs/K^m ^ Ao + Bo (M^/Mgoo)"", in the 
range 0.1r2oo — ^500- We give mean and rms of the best 
fit values in Table [TJ Note that the power-law indices 
for the observed entropy-gas mass relation are shallower 
than expected from the theoretical expectation. Interest- 
ingly this index (logarithmic slope) does not differ much 
in the whole cluster sample. The values of Ao and Bo 
show that cool-core clusters are entropy deficient. 

4. EFFECTS OF ONLY PREHEATING AND COOLING 

iVoit et ahl ()2002D discussed three types of modifica- 
tions to the initial entropy profile, namely: (i) a trun- 
cation in the entropy profile owing to removal of gas, 
approximating the effect of gas cooling and dropping out 
of the ICM, (ii) a shift in the profile, mimicking the ef- 
fect of preheating, and (iii) lowering the entropy profile 
due to radiative cooling. Assuming a form of the cooling 
function of the type A oc T~^ " for group temperatures 

(T < 2 keV), it was shown that K^^^ across the cluster 

.5 3/2 
is reduced by an amoun t ^Kc , where Kc is a critical 

entropy. iJohnson et al.l (|2009( ) suggested a combination 
of the effects of preheating and cooling, expressed as 
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th 



K 



shift 



l3/2 _ £&'3/2 
2 '^ 



(5) 



Kc « 81keVcm2[T/lkeV]2/3[f/14Gyr]2/3, (their eqn 
14) describes the cooling. They calculated the constant 
preheating shift K shift by evaluating eqn [5] at their out- 
ermost radial point. 




Fig. 2. — The ratio of the theoretical to observed temperature 
as a function of Mg/M^oo is plotted for our sample of 25 clusters. 



For the 28 n earby galaxy g roups from the XMM- 
Newton survey, IJohnson et al.l ()2009D found that while 
this 'preheating + cooling' model matches the observa- 
tions better than a simple shift/truncation, it still fell 
short of being a reasonable representation of the observed 
profiles. 

We have fitted models of the form of eqn [5] to our 
sample, where we have used the entire profile Kth{Mg) 
and Kobs{Mg), rather than just one radial point for the 
fit. We attempted three different types of fits for each 
cluster, described below: (i) Kc is evaluated using the 
full radial temperature profile instead of mean temper- 
ature; (ii) A fit using the constant Tgp for each clus- 
ter; (iii) We assume that a fraction of the gas mass is 
lost from the ICM, and try two fits with varying frac- 
tions of the total gas mass, / = 0.8, 0.9. The tempera- 
ture T = Tsp, and the expression used for the fitting is: 

KobsiJMg) = [{Kth{Mg) + Ksh^ftf'^ - (3/2)i^5/']' ^ 

We find that the number of clusters for which none of 
fits are good far exceeds the clusters for which any of fits 
can be called reasonable (reduced x^ < 2). The lack of a 
good fit to the preheating+cooling model in most of the 
clusters in the sample suggests that a major component 
of entropy enhancement occurs beyond simple preheating 
and radiative cooling. 

5. FRACTIONAL ENTROPY DEVIATION AND ENERGY 
INPUT 

In order to determine the amount of energy deposi- 
tion associated with the entropy enhancement, we use 
the quantity Tobs^K/Kobs, where A/C = Kobs - Kth- 
This is because the thermodynamic entropy of an ideal 
gas S is related to K, as S = Const. x\nK and the 
change in energy per unit mass dO = Td S oc TAK/K 
(see also eqn 3 of iFinoguenov et aLl (120081)') 



In an isochoric process, (see also (jLlovd-davies et all 
[2000)), 



AQ 



AKn. 



2/3 



kTob. AK 



(7 - 1) (7 - 1) Kob 



(6) 
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Fig. 3. — The profiles of energy deposition per unit particle 
plotted against M9/M500, after scaling them by Tsp for different 
clusters. Data points inside the core radii are shown in dotted lines. 
The color scheme is same as that used in fig[T] The mean profile 
and 1-fT scatter, outside the core radii, is shown with thick black 
line. 

In an isobaric process, however, (for Tf/Ti — /3) 



AQ = 



AKn 



2/3 
/ 



(/?' 



5/3 



1) 



(l-i)Mm,/32/3(/3-l) 
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The ratio of the changes in energy for a given fractional 
change AK/Kobs and Tots, is given by. 



AQ 
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AQ 
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isocnortc 



/32/3(^_l) 



(8) 



For a value of /3 = 2, the ratio is 1.14. This implies that 
if the observed temperature To^siMg) deviates from the 
theoretically calculated value Tth{Mg) by a factor < 2, 
then the two above mentioned estimates of energy input 
per unit mass differ by only a factor of 1.14. Figure [2] 

shows the ratio ^"'V,/ ^ for clusters in the sample. One 

can see that the two temperature profiles vary within a 
factor of ~ 2 ; hence, we choose the expression for the 
isochoric process in our estimates. 

We first estimate the energy per particle, AE{Mg) = 
^Tobs^^, for each cluster. Fig [3] shows the profiles 
for AE/Tsp , the ratio of the non-gravitational energy 
injection to the gravitational potential of the clusters in 
three temperature bins. While the detailed profiles differ 
from cluster to cluster, AE/Tgp generally has a decreas- 
ing profile with a similar trend. Albeit with a large scat- 
ter within each temperature group, we find that the mean 
value of AE/Tsp is the highest for the lowest temperature 
group and vice versa. For all clusters, non-gravitational 
energy is already comparable to gravitational energy at 
core radii. Moreover, the profiles decrease by a factor 
of - 50% for Tsp < 3.5 keV clusters and - 75% for 
Tsp > 3.5 keV clusters. Thus, on average, the profiles 
for higher masses decrease faster than those for lower 



masses. Our calculations, as mentioned earlier, are valid 
outside the core 0.1r2oo- 

We determine the mean profile after averaging over all 
the AE profile fits for the individual clusters^ The mean 
profile with the 1 — a scatter is shown in figE]. The mean 
profile decreases by roughly a factor of 4 between 0.1r2oo 
and rsoo. 

The total amount of energy deposited, for the whole 
cluster is , 



^non — grav 



T^^^dMg, (9) 

(7 - l)fimp Kobs 



for Mg/M^oo between the hmits 0.1r2oo < r < r^oQ. We 
use the fits obtained above to calculate the integral and 
the results are shown in Figure ID Clearly Enon-grav 
is proportional to the cluster mass. A fit results in the 
following scaling relations for the whole sample: 



E^, 



non — grav 

10^1 keV 



(-0.414±0.41)- 



-(0.2±0.17)(^)i-6^±°-4^ 
keV 



(10) 

If the cool-core clusters are omitted, then one obtains a 
slope of 1.53 ±0.64, a constant term of —0.48 ±0.65 and 
first coefficient as 0.24 ± 0.3. 

Dividing the energy by the total number of particles in 
the ICM, we estimate the mean energy to be 2.74 ± 0.87 
keV per particle. 
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Fig. 4. — The total energy injection between 0.1r200 — ^500 is 
plotted against cluster mass. Red points are for non cool-core 
clusters and blue points show cool-core clusters. The best fit to 
entire sample is shown in solid black line. The black dashed line 
shows best fit to the non cool-core clusters. 



6. DISCUSSIONS 

We can compare the value of ~ 2.74 keV per parti- 
cle as a measure of the energy input corresponding to 
the observed entropy enhancement, with those inferred 
from other considerations. We note that the earlier ob- 
servations of entropy enhancement in the ICM, mostly 
inferred from the deviations of several cluster scaling re- 
lations, were interpreted in terms of an entropy floor. It 
was thought that elevating the ICM gas to an entropy 

3 To this end, we use the fit A_B = C + D (Mg/Msoo)* 
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floor of a few hundred keV cm^ would explain the ob- 
servations. This entropy floor can be associated with 
an amount of energy if the density of the gas at the 
epoch of energy input i s know n. For example, simula- 
tions by IBorgani et al.l (|2001D showed that an entropy 
floor of ~ 50 keV cm^ was adequate to explain the ob- 
servations. They found this corresponds to ^ 1 keV per 
particle if the heating was assumed to have taken place 
at z '-^ 3. 

Interestingly, iRovchowdhurv et al.l ()2004[ ) showed for 
their model of AGN feedback from black holes that X- 
ray observations could be explained with an energy input 
proportional to cluster mass. In their model of AGN 
feedback through buoyant bubbles of relativistic plasma, 
which deposit energy into the ICM through pdV work, 
convection and thermal conduction, this proportionality 
implied a linear relation between the black hole mass of 
the central AGN and the cluster mass. We find that in 
order to explain the correlation in Figure IH we need the 
black hole mass Mbh ^ 2 x 10~^M5oo ?yo.2, where the 
energy available from the AGN is characterized by an 
efficiency rj = 0.2^ 

Simulations bv lGaspari et al.l ()2012D also show that the 
energy deposition is centrally peaked within the core of 
the cluster. We note however that our results pertain 
only to regions outside the core radius. 

We have estimated the energy input corresponding 
to the entropy enhancement differently from previous 
works: Firstly, we have not used any cluster scaling rela- 
tions which depend on the average properties of the ICM. 
We have also used the distribution of the X-ray entropy 
(K) with gas mass, since entropy per unit mass (5*) is a 
Lagrangian quantity. Furthermore, instead of determin- 
ing an entropy floor and then estimating an amount of 
energy assuming a certain density, we have estimated the 
energy input from first principles. 

Our result implies that the effect of energy deposition 
in low and high temperature clusters is remarkably simi- 
lar although the entropy enhancement processes may dif- 
fer substantially in these systems. Secondly, the gas mass 
profiles of energy deposition per particle show that the 
processes responsible for entropy enhancement in clus- 
ters affect the gas in the central regions more than in 
the outer regions. Moreover, profiles for higher masses 
decreases faster than those for lower masses. The ba- 
sic similarity (modulo some dependence on the cluster 
masses) in the profiles can provide a test for future sim- 
ulations. 

Taken together, our results indicate that the energy 
associated with entropy enhancement is proportional to 
cluster mass. Furthermore, their effect in all clusters is 
centrally peaked. This suggests an energy source which 
must satisfy both requirements simultaneously. As men- 
tioned earli er, AGN feedback models s atisfy the first re- 
quirement (jRovchowdhurv et al.ll2004 ). It is also plau- 
sible that the effect of the feedback is more pronounced 
in the inner regions, driving most of its gas outside th e 
inner region (jMcCarthv et al.ll2010l : lGasDari et al.ll2012[ ). 

7. SUMMARY 

We have looked at the fractional entropy enhancement 
in the ICM for a sample of REXCESS clusters by compar- 
ing the observed entropy profiles to those expected from 
gravitational collapse only. We first show that this en- 



tropy excess cannot be explained by only preheating plus 
cooling models of entropy enhancement. Since, this en- 
tropy excess must be sourced from non-gravitational pro- 
cesses, we connect this excess to any non-gravitational 
energy deposition in the ICM. We report, to our knowl- 
edge, the first energy deposition profiles in a large sample 
of clusters and also estimate the total non-gravitational 
energy that has been dumped into the ICM. We find that 
this excess energy is proportional to cluster temperature 
(and hence cluster mass). We show that the entropy en- 
hancement process in the ICM is centrally peaked and 
is relatively larger in low temperature clusters than in 
high temperature clusters. Our results support models 
of entropy enhancement through AGN feedback. 
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